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Abstract 

The interface between the insulating oxides LaAlOs and SrTiOs exhibits a superconducting 
two-dimensional electron system that can be modulated by a gate voltage. While gating of the 
conductivity has been probed extensively and gating of the superconducting critical temperature 
has been demonstrated, the question whether, and if so how, the gate tunes the superfluid density 
and superconducting order parameter is unanswered. We present local magnetic susceptibility, 
related to the superfluid density, as a function of temperature, gate voltage and location. We show 
that the temperature dependence of the superfluid density at different gate voltages collapse to a 
single curve characteristic of a full superconducting gap. Further, we show that the dipole moments 
observed in this system are not modulated by the gate voltage. 

PACS numbers: 74.78.Fk, 71.70.Ej, 74.25.N-, 75.20.-g 
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Electric field control of conducting channels has allowed great innovation in traditional 
semiconductor devices l| . Now heterointerfaces in a new class of materials, the complex ox- 
ides, have generated significant interest because of their gate tunable properties. Specifically, 
the conducting interface formed between the band insulators lanthanum aluminate and TiOo 
terminated 100 strontium titanate (LAO/STO) [2] exhibits many fascinating properties js] 
suggesting that an electronic reconstruction triggered by the polar/non-polar interface plays 
an important role in the inducing the conductivity in the STO Q . At low temperatures this 
interface displays two-dimensional superconductivity 5^. 



constant of STO at low temperatures 



Additionally, the high dielectric 



6| makes applying an electric field with a back gate 
especially effective to tune the properties of this superconducting state. 

Caviglia etal. showed that with increasing gate voltage, Vg, the superconducting critical 
temperature, Tc, displayed a dome structure and concurrently the normal state resistance 
monotonically decreased ^. Later work showed that the electron mobility and carrier 
density both increased continuously with Vg, with the former dominating the Vg dependence 



of the conductivity jsj. The evolution of a non-linearity in the Hall resistivity as a function of 



9| has been interpreted b y J oshua et al. as evidence of electrons populating conduction 



10|, implying that the ratio of high and low mobility electrons 



bands with different mobilities 
may be tuned by gating. 

Notably, the interface breaks spatial inversion symmetry, opening the possibility for spin 
orbit coupling to impact the electronic properties of the interface gas. Two groups reported 
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tuning of the Rashba spin orbit coupling (RSOC) inferred from magnetoresistance 
and measurements of the in-plane critical fields |9|. They found opposite dependencies for 
tuning the strength of the spin orbit coupling with Vg, making the impact of Vg on the spin 
orbit coupling unclear, possibly suggesting a peak in the spin orbit coupling. 

Moreover, the discovery of magnetic patches coexistent with superconductivity 
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-ll5| 



and the presence of RSOC originating from the noncentrosymmetric nature of the interface 
have raised the possibility of an unconventional superconducting pairing mechanism or order 



parameter 



16l4l8l|. However, all previous measurements studying how gating effects the 



properties of the interface used electronic transport, which gives limited information about 
the superconducting state. In this Letter, we use local magnetic susceptibility to make the 
first direct measurements of the superfiuid density in LAO/STO and address the question 
of how the superconducting state evolves with Vg. 



Measurements were made on a sample with five unit cells of LAO grown at 800° C and 
1.3 X lO^^mbar oxygen partial pressure on a Ti02 terminated STO substrate. The growth 
was followed by a high pressure oxygen anneal, 600° C in 0.4 bar. The sample was silver 
epoxied to a piece of copper tape, which served as a back gate. Vg was applied between the 
copper tape and the interface, which was contacted by aluminum wirebonds. Magnetiza- 
tion and susceptibility measurements were made using a scanning SQUID (Superconducting 
Quantum Interference Device) [isj] , with a 3 /im diameter pick-up loop and a concentric field 
coil for applying a local AC magnetic field. The pick-up loop is sensitive to both the DC 
static flux and the AC flux resulting from diamagnetic screening currents cancelling the field 
from the field coil. This setup enables simultaneous measurements of ferromagnetism and 



superconductivity in the sample 
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A superconductor will generate screening currents to screen an applied field. The currents 
extend into a bulk superconductor by the penetration depth, A. The temperature dependence 
of A is a probe of the superconducting state. For a thin superconductor of thickness d, the 



screening distance is given by the Pearl length A = 2\^/d ^0|. Using a model by Kogan 
{l| , we extract A from measurements of the screening currents as a function of the distance 
between the sensor and the sample. A is related to the superfiuid density, Ug = 2m*//ioe^A, 
where e is the elementary charge, /xq the permeability of free space, and m* = lAQrrie the 
effective electron mass measured by 22| from Shubnikov de Haas on LAO/STO interfaces. 
We repeat these measurements at multiple temperatures and gate voltages to map out the 
superconducting state. Fig. 1. We define Tc as the temperature at which the diamagnetic 
screening drops below our noise level of 0.01 $o/A, corresponding to a minimum detectable 
of 4 — 14 X 10^°cm~^. The statistical errors were smaller than the systematic errors. 



outlined in gray in Fig. 2a, from imprecise knowledge of our measurement geometry 23 |. 
The systematic errors are fixed for a single cooldown and represent an overall scaling of 
which would be the same for every measurement. 

Tc vs Vg (Fig. la) has a maximum Tc = 240 mK. In the range of applied Vg super- 
conductivity can only be eliminated on the underdoped side of the dome, and Ug grows 
monotonically with Vg, with Ug = 3.0 x 10^^ cm~^ at the largest Vg. (Fig. lb) The carrier 
density and mobility were measured in a separate cooldown with no backgate. At 2 K the 
mobility was 1.02 x lO^cm^/Vs and the density was 2.05 x 10^^ cm^^, ten times larger than 
the largest Ug we observed. 
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FIG. 1. a) The critical temperature as a function of gate voltage forms a dome. The dashed line 
represents our lowest measurement temperature, b) The superfluid density at our lowest tempera- 
ture as a function of gate voltage. The superfluid density increases monotonically throughout the 
dome. The color scale represents gate voltage and is repeated in Fig. 2. 

A small ratio of the superfluid density to the normal density is expected in the dirty 
limit, in which the elastic scattering time, r, much shorter than the superconducting gap, 
Aq (h/r ^ Ao). h is reduced Plank's constant. Above the normal density of electrons n is 
given by the optical sum rule n oc ai{uj)du, where o"i is the real part of the conductivity 
and u the frequency. For a metal a is sharply peaked near zero frequency, so scattering 
moves spectral weight to higher frequencies. Below Tc, a gap opens at u; = 2Ao/h and the 
spectral weight within that gap collapses to a delta function at the origin whose amplitude 



is proportional to Ug 



24j |. Therefore in the dirty limit, only a fraction of carriers enter the 



superconducting state, n^/n = 2Ao/(^/r). Using the gate tuned mobilities reported by Bell 
etal., 100 — lOOOcm^/Vs [sl, we expect the ratio Us/n to be 0.01 — 0.1, consistent with our 
measured n^. 

We now look at the temperature dependence of the superfluid density. Fig. 2b plots Ug vs. 
T for all Vg across the dome. Strikingly, when normalizing the curves they collapse (See Fig. 
2c), showing that within our experimental errors there is no change in the superconducting 
gap structure with electrostatic doping. Furthermore, the collapse is reproducible over 
multiple positions, sweeps of Vg, and samples [25|. 

The temperature dependence of the superfluid density is a direct probe of the supercon- 
ducting order parameter. It can be used to distinguish BCS superconductors from unconven- 
tional superconductivity. We fit^the normalized curves to a phenomenological BCS model 
with two parameters A and a [2,]. A scales the superconducting gap Aq = A/c^Tc. a is a 




FIG. 2. a) Superfluid density vs. temperature for Vg = llOy, the peak of the superconducting 
dome. The gray area shows systematic error, b) Superfluid density vs. temperature for every gate 
voltage. The colors represent the same Vg from Fig. 1. c) Normalized curves from b). The gray 
line shows the temperature dependence of a weakly interacting clean BCS s-wave superconductor 
(A = 1.76 and a = 1). The black dashed line is a fit to the data (A = 2.2 and a = 1.4). 



shape parameter that determines how rapidly the gap opens below Tc, Us oc 1 — {T/Tc 
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27l |. A = 1.76 and a = 1 for an clean s-wave BCS superconductor with weak coupling 
2|, plotted as the gray line in Fig. 2c. The fit to our data gives A = 2.2 and a = 1.4. 
This is consistent with a BCS description with increased coupling or disorder. Both will 



theoretically increase the gap and the a parameter 281], shifting the curve up and to the 
right. 

The flattening at low temperature indicates fully gapped behavior with a gap that is 
larger than BCS weak-coupling s-wave. Our lowest measurement temperature is 1/6 of 
j-max^ and Hs remains fiat (within 3%) up to 35% of T^. A full gap indicates the absence of 
low energy quasiparticle excitations, ruling out order parameters with nodes in the Fermi 
surface. Furthermore, the steep rise of Ug near Tc and the absence of a kink in the functional 



form rule out most weak coupling two band models [29|, because a second smaller gap will 



slow the onset of superconductivity near T^. Two gaps of similar size, both larger than the 
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FIG. 3. Critical temperature vs. the superfluid density at lowest temperature (T ~ 40 mK). 
The red points are the data from Figure 1 and the gray dots represent additional data sets. The 
dotted line is the theoretical phase fluctuation temperature from ref [30(], which may be limiting 



the critical temperature on the underdoped side of the dome. The bimodal distribution on the 
overdoped side is due to inhomogeneity that locally suppresses Ug in different regions of the sample 
while the Tc remains the same. See also Fig. 4. 



BCS gap or a dominant single large gap with second smaller amplitude gap, could reproduce 
the data. 

The low Ug in the underdoped region may result in suppression of Tc by thermal phase 
fluctuations. Such fluctuations would result in a linear temperature dependence of in 
the underdoped region. Following reference 301], we calculate a phase ordering temperature, 
rpmax _ j\fi'^fi^(^Qy 4fi2* ^ where A = 0.9 in two dimensional systems. Fig. 3 shows Tc vs 
ns(40mK), additionally TJ""^ is plotted as a linear function of n^: the line does not suggest 
a fit to our data. We have insufficient data at the lowest superfluid densities to make any 
statement about the functional form of Tc(ns) in the region where phase fluctuations may 
be limiting T^. Nevertheless, the proximity of the phase ordering line to the underdoped 
data suggests that phase fluctuations may drive the abrupt decrease of Tc. 

Given the 2D nature of the superconducting system we expect a BKT transition, where 
unbinding of vortex anti- vortex pairs suppresses superconductivity and results in a discontin- 
uous jump in Hg near T^. The jump should occur at flnite superfluid density = 2m*Tc/ nh'^ 



3l| . For the maximum Tc = 240 mK a BKT transition should occur at 5 x 10^*^ cm ^, which 



is too close to our measurement threshold to establish a BKT jump in our Ug vs. T curves. 
Are our observations consistent with a simple s— wave order parameter from doped STO 



32l | or a two gap mixed state induced by symmetry breaking at the interface? Rashba spin 
orbit coupling (RSOC), induced by the structural inversion asymmetry, is expected to lift 
the spin degeneracy and split the energy bands [33]. Additionally, RSOC breaks parity 
and consequently mixes singlet and triplet states resultin g in an s— wave component 



mixed with a triplet induced d- vector d(k) = xky — yk^ 18|, |3J]. Mixing results in two 
gaps, A = As ± \dk\, whose magnitudes depend on the weights of the singlet and triplet 
components. Varying the relative weights changes the density of states, but always results 
in two fully gapped Fermi surfaces except for the special case where the s-wave singlet and 
triplet gaps are the same and accidental line nodes form on one band [isl l. 

Other reports [o], [lli have demonstrated significant tuning of the strength of RSOC with 
Vg. An open question, of particular importance to testing this two gap picture, is how do the 
weight of the two components change with Vg. Our results, showing a consistent functional 
form for vs. T across all Vg, show that the superconducting gap structure does not change 
with Vg consequently the relative gap weights do not change with Vg. The effect of RSOC on 
the band structure may depend on the chemical potential which is also tuned by the gate. 
Therefore the insensitivity of superconductivity to Vg cannot completely rule out a RSOC 
induced two gap scenario. Yet, our second observation of the fast opening of the gap near 
Tc and the compatibility of the data with a single gap BCS model limits two gap models. 
Both gaps must be larger than the BCS s-wave gap to capture both the fast rise and flat 



low temperature dependence of the data [35 1. 



Finally, disorder may play a role in washing out the triplet component. As stated above, 
the LAO/STO system is a dirty superconductor, with h/r » A. Disorder averaging 
has very little impact on the isotropic s-wave component but may eliminate the triplet 
component. 

In short, our data is most consistent with a single gap. We cannot rule out the presence 
of two gaps, but our observations limit their size and Vg dependence. 

Our scanning SQUID system allows two dimensional mapping of superconductivity and 
magnetism at different Vg. Fig. 4 shows simultaneously imaged susceptometry and magne- 
tometry scans of the same region at 80 mK for four different Vg . The inhomogeneity in the 
diamagnetic screening is very large in the underdoped region [Vg = — lOV) and re-enters 
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FIG. 4. Susceptometry (left) and magnetometry (right) at 80mK at different gate voltages. 
(Inset) Reproduction of the Tc dome from FIG 1 showing the relative location of Vg in each 
panel, a-b) The sample is no longer superconducting and has a paramagnetic response. Individual 
ferromagnetic dipoles are also visible in the paramagnetic image, c-d) Superconductivity appears 
and the landscape is relatively inhomogeneous. e-f) Peak of the superconducting dome, most 
inhomogeneity disappears, g-h) Excess inhomogeneity returns on the overdoped side of the dome. 
The ferromagnetic patches do not change with Vg and remain when superconductivity is gone. 

the image in the overdoped region {Vg = 390 V). The least inhomogeneity is observed at 
optimal doping, although it does not disappear. In contrast the ferromagnetic patches are 
insensitive to Vg with a constant magnitude and orientation for all Vg. This behaviour was 
also observed on 15, 10 and 3.3 uc samples, showing the electron density that is modified 
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by Vg does not appear to influence the ferromagnetism. 

In conclusion, we presented the flrst measurements of the superfluid density as a function 
of temperature at multiple gate voltages throughout the superconducting dome in LAO/STO 
heterostructures. The temperature dependence of is well described by a fully gapped 
BCS model. Moreover, the normalized Ug vs. T curves collapse to a single functional form 
indicating there is no change in the gap structure with Vg. Although we cannot rule out a 
two gap mixed singlet /triplet model, the insensitivity of the superconducting state to Vg and 
the large slope near Tc limit two gap scenarios. Speciflcally, both gaps must be larger than 
the BCS s-wave gap and their relative size cannot change throughout the dome. A future 
experiment to distinguish between these two scenarios may be to gate the superconductivity 
in the presence of an in-plane fleld, which can change the relative magnitude of triplet and 
singlet gaps. Alternatively, samples in the clean limit may reveal a clearer two gap structure. 
Additionally, we found that the magnitude and orientation of the ferromagnetic patches that 
coexist with superconductivity are unchanged by Vg, while at the same time Ug goes from 
zero to 3.0 x 10^^ cm~^. This shows the population of electrons that is modifled by the gate 
is separate from the electrons that contribute to the ferromagnetic order. 
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SOM: Gate tuned superfluid density measurements of supercon- 
ducting LaAlOg/SrTiOa 



DISCUSSION OF SYSTEMATIC ERRORS 

The accuracy of our superfluid density measurement is dominated by systematic errors 
wfiicli arise from insufficient knowledge of the physical parameters of our SQUID sensor and 
piezoelectric scanner. Our measurement of the superfluid density, n^, relies on extracting 
the Pearl length, A, from fits to approach curves. 

2m* . . 

An approach curve measures the diamagnetic susceptibility as a function of the sensor height 
above the sample. Our SQUID sensor consists of a pair of concentric current carrying wires 
called the field coil and pick-up loop. 

shields 





lOum 



FIG. 5. a) Actual layout of the SQUID field coil, pick-up loop and shields, b) Approximations to 
the actual layout used by Kogan 

We follow a model developed by Kogan which treats the SQUID's field coil as a circular 
current loop of radius, a When the loop is brought near a superconducting thin film, 
the Meissner response of the film detected by the pick-up loop can be expressed as 

= /iovrapy^ dk^-^e-^"'J,{ka)J,ikp), (2) 

where p is the pick-up loop radius, and /io is the magnetic constant. This integral gives a 
value for the diamagnetic susceptibility, $, at a height h above the sample. Six physical 
parameters enter equation ([2]): the radius of the pick-up loop p, the radius of the field coil 
a, the piezo calibration from volts to microns Vc, the distance between the pick-up loop and 
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the sample when the SQUID makes contact ho, the offset of the susceptibility far from the 
sample ^off, and a background slope m. We convert the voltage applied to the z-bender, 
T4, to a height h — T^K + ^o- The susceptibihty seen by the SQUID is 

^SQ^^ + ^off + mh. (3) 

Consequently, our fits for A depend on the accuracy of our knowledge of the other parameters. 

We start with estimates of p and a. We can make accurate measurements of the two radii 
using an optical microscope; however these wire loops have a finite width and leads that 
deform their magnetic response with respect to the perfectly circular loops in Kogan's model. 
Using numerical methods, we calculate the source field using the measured dimensions of 
our non-ideal coils. We find the non-ideal nature of the field coil and pick-up loop results in 
a 15% error on the fitted value of A. This error works out to few hundred microns on our 
shortest A fits. 

We now address the errors associated with our bender constant Vc and height offset ho. 
We don't have accurate calibrations for these parameters, but we do know that these values 
should be the same for every touchdown curve. We fit hundreds of approach curves using 
A, Vc, and ho as free parameters, and assembled histograms of the fitted Vc and ho values. 
Prom the histograms we were able to extract a best value and variance, a. We then use the 
error propagation equation to relate the variances in Vc and ho to an error in A. 

^ < (t)' ^ < ■ + (t) {IQ ^ ■ 

The propagation equation yielded an error of about 1 mm on our shortest A fits. This is a 
systematic error and is the same for every touchdown curve in the cooldown. It may change 
the overall calibration for Us, but it will not change the trends in Us vs or vs T. 

We added the systematic errors from the sensor coils, bender calibration and height offset. 
The total systematic error is show as the gray outline shown in Fig. 2a of the main text. 
The error from the bender and offset dominates the error from the non-ideal nature of the 
pick-up loop and field coil. 

DISCUSSION OF PHENOMENOLOGICAL BCS FITS 

We compare our normalized plots of superfiuid density vs. temperature to a phenomeno- 
logical BCS model with an isotropic s-wave superconducting gap. The normalized superfiuid 
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density, ns/ns{T = 0), was given by Prozorov and Giannetta 



r = 1 - — / cosh"^ ^ ^ de, (5) 

ns{T = 0) 2T io \^ 2T y ' ^ ^ 

where T is the temperature and A(T) is the superconducting gap function. The gap can be 
written Isl as 



A„(r) = A„(0)ta„h^^^a(|-l)j. (6) 

Ao(0) is the zero temperature energy gap and a is a shape parameter which determines how 
fast the gap opens. Near the critical temperature the superfluid density can be approximated 
as ns = 1 — {T /TcY°-. For an isotropic s-wave gap Ao(0) = l.TG/c^Tc and a = 1. We use 
equations (|5]) and ([6]) to fit our data and find Ao(0) = 2.2A;bTc and a = 1.4. This is the 
dashed hne plotted with the data in Fig. 2c of the main text. 



DISCUSSION OF TWO GAPS IN BCS 



We can use equations ([5]) and ([6]) to generate a phenomenological two-gap expression [4 1 . 

n,(T) = pn,i{T) + (1 - p)n,2{T) (7) 

Fig. 6 shows plots of the superfluid density for two gaps of equal weight {p = .5) with 
different physical parameters. The only combination that can support a superfluid density 
function that rises faster than BCS near Tc has two gaps that are larger than the BCS gap. 
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FIG. 6. Comparison of a two-gap superfluid density Ug with a single gap BCS superfluid density 
^sBCS- In ^-ll three plots n^i = pusi and 71^2 = pns2 with a) Plots of two gaps with Ai = A2 = 1.76 
and ai = a2 = 1 but two different critical temperatures, b) Plots of two gaps with Ai = 2.2, 
A2 = 1.1, and oi = 02 = 1. c) Plots of two gaps with with Ai = 3, A2 = 2, ai = 1.8 and 02 = 1. 
Only in c) where both gaps are larger than Abcs = 1.76 can we generate a total superfluid density 
that opens faster than BCS. 
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